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t;;^- , The two-dimensional uniformly driven diffusive system is simulated with a multi-spin 
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0^ ! coding technique. The nonequilibrium phase transition is analyzed with anisotropic 
^ , 

Q ' finite-size scaling, both at the critical point and off the critical point. The field- 

D ; 

Q ' theoretical values of critical exponents fit the data well. 
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'nI" ! 1. Introduction 
5 ■ 

^ • Driven diffusive systems are a class of models which exhibit nonequilibrium phase 

a : 

I . transitions. The models are defined by some local rules similar to Kawasaki dy- 

^ : 

C , namics. The steady states of the dynamical evolution have been studied extensively 

' during past decade. A central issue is whether the concept of universality of critical 
> . 

I phenomena can be applied in nonequilibrium cases. 
c3 ' The standard driven diffusive model of half-filled charged lattice gas was proposed 

as a model for ionic solution in an electric field. A continuum version, based on 
symmetries and conservation laws, was solved in a field-theoretic framework. It is 
quite remarkable that exact critical exponents are obtained for dimensions from two to 
five. In particular, the set of critical exponents in two dimensions is 

1 3 1 

These exponents have similar meanings as in equilibrium second-order phase transi- 
tions. They are thought to be universal within a class. No independent methods have 
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been applied to the problem to check out the validity of the claim except computer 
simulations. C'''^) Anisotropic finite-size scaling studies by Leung(^'-^°) appear to have 
settled a dispute between theory and early computer simulations. But a recent work^^^ 
casts doubt on Leung's conclusion. We present here high quality data to support Le- 
ung's results. 

2. Model and Simulation 

The system consists of a square lattice of L x M sites. The driven field is in x 
direction. A configuration is specified by a set of Ising spins, ax^y = ±1, with zero 
total magnetization. The state evolves according to the following prescription. A 
bond is chosen with equal probability in orientations and in locations. If the bond is 
parallel to the driven field and the adjoining spins are distinct, the spins are changed to 
(^x,y = ~1 cind (Tx+i mod L,y = +1- This Corresponds to an infinitely strong driven field. 
If the bond is perpendicular to the field, we swap the spin values with a probability 
minjl, exp(— 5£'//cbT)}, where SE is the energy increment due to the change, assuming 
the usual nearest neighbor ferromagnetic interaction with coupling constant J and 
periodic boundary conditions. One Monte Carlo step is defined as L x M such basic 
steps. 

Simulation near critical point is not easy because of critical slowing down. The 
situation for this model is more severe by the conservative nature of the dynamics, 
leading to relaxation time r oc M^. Thus, an efficient implementation is crucial to 
obtain good statistics. With a multi-spin coding method, '■''^''^^ thirty two systems are 
simulated simultaneously. This gives us at least a factor of 10 speed up. A slight 
penalty of the multi-spin coding algorithm is that the temperature can not be set 
exactly to the desired value. But it can be well under control with a large random bit 
table. To achieve an accuracy of five significant figures, we took a random bit table of 
2-'^^ entries. 

The program runs at 0.3 //sec per spin exchange on an SGI Iris Indigo Workstation. 
Computations were done on a cluster of over fifty workstations (SGI Iris Indigo, HP 
9000/700 model 712, and DEC 5000) in two months. The lengths of the runs are lO'^ 
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to 10^ Monte Carlo steps. These are orders of magnitude longer than previous studies. 
Measurements are performed at an interval of 10 Monte Carlo steps. 

We use the order parameter introduced by Binder and Wang^^^^ and modified by 
Leung. Let's define 

L-lM-l 

/ TT \ 
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The normalization is such that (/) = 1 for a strip geometry (the configuration in the limit 
T ^ 0). The following quantities are calculated, (1) the order parameter ^' = (0), (2) 
the "susceptibility," or fiuctuation of the order parameter. 
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and the susceptibility above the critical temperature, 
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and (3) the fourth-order cumulant. 



g = 2- 
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Note that g goes from to 1 as temperature T goes from oo to 0. We'll measure 
temperature in units of the two-dimensional Ising critical temperature (2.269J/A;b). 

3. Determination of the Critical Temperature 

Estimating the critical exponents depends crucially on an accurate value of critical 
temperature. In the literature different authors have given incompatible values, = 
1.355 ± 0.003 by Valles and Marro,(^) = 1.30 ± 0.01 by Achahbar and Marro,^^) and 
Leung's result*^^-' of = 1.418 ± 0.005. The discrepancies are the manifestation of the 
difficulties of simulating the system and of interpreting the complicated finite-size data. 

We determine using two independent methods. First, we look at the peak of the 
susceptibility x- The method is not particularly accurate, but we get a general picture 
of the finite-size critical temperature Tc{L,M). For very elongated shapes, Tc{L,M) 
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decreases towards zero. For a fixed value M or L, Tc{L, M) reaches its maximum at 
about S = V-/^ /M Ri 0.2. The systems with ratio S 0.2 notice the finiteness of the 
sizes in two directions at roughly the same temperature. Tc(L, M) increases as system 
size increases with fixed S. Figure 1 is a plot of the peak locations Tc(L, M) versus 
L~^/^, for some fixed values S as well as for the square geometry L = M. According 
to the usual finite-size scaling assumption, we should expect 

Te(L,M) = Tc(oo,oo) + F(5)L-^/"ii. (6) 

Clearly, the data do not follow this equation very well, presumably due to corrections 
to scaling. But linear extrapolations should give us lower bounds for the infinite system 
transition temperature Tc. We quote the following result, 

Tc = 1.410 ± 0.006. (7) 

It agrees with Leung's result*^^-* within errors. For the square systems, the convergence 
to the critical temperature is extremely slow. This may explain why the previous 
calculations on square systems all gave a lower Tc. If anisotropy is a dominant feature, 
we should not expect peak locations to scale simply as L~^l^ for square systems. 

The peaks were not located with great precise because the simulations were carried 
out at discrete points spaced at AT = 0.01. The second method exploits the scaling 
properties of the fourth-order cumulant. From finite-size scaling theory, we have 

^(T, L, M) = ~g{l}ln (T - Tc)/Tc, S) . (8) 

If scaling were exactly obeyed, different curves (with fixed S) should intersect at exactly 
the same value Tc. Therefore, there is no need to extrapolate to infinite size. In practice 
there are unknown corrections to scaling. A better way is to consider the overall scaling, 
Eq. 8. The value of Tg for each S can be determined more precisely. However, there 
are weak size and S dependence. Nevertheless, we found that the values all fall into 
the interval 1.395 to 1.410. Figure 2 is one of the scaling plots with S = 2~^/^. 

Even with much greater computational effort, the accuracy of Tg has not been 
improved. Each set of data or each method may give more accurate value, but diflFerent 
sets of data or methods give slightly different values. 
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4. Anisotropic Finite-Size Scaling at Critical Point 

The finite-size scaling theory of isotropic systems can be generahzed to anisotropic 
systems. The idea is that the anisotropic system has the same scahng form as the 
isotropic one if we fix the ratio S = L^-^/^ii /M. In anisotropic finite-size scahng, S enters 
as an independent variable, in addition to the usual scaling variable L^^'^w (T — Tc)/Tc. 
particularly, at T = Tc, scaling functions depend on S. 

In applying the finite-size scaling theory, we could simulate a very large system and 
study the size effect of smaller subsystems. This may seem computationally effective. 
But there are two problems associated with it: we have the annoying finite-size effect of 
the very large system when the subsystem sizes are comparable to it; and we may not be 
able to equilibrate the very large system very well. So, we adopt the more conventional 
finite-size scaling analysis — working on fully finite size of dimension L x M. 

The exponent ratio i^\\/i^± is one of the most important number in an anisotropic 
finite-size analysis. Theoretical result is often assumed. '^^'■'^^^ We have attempted to 
determine it from our data. The strip geometry with — > or — > oo and periodic 
boundary conditions has a simpler scaling behavior at the critical temperature, '^^^^ 

x{Tc) (xM''/"^, *(Tc) oc L>M^ii/^^, (9) 

x(Te) (X L^/^ii , *(Te) (X L^/^^ii-^/^M-i/^ M>L'^^/^ii. (10) 

Figure 3 shows the long strip limiting behavior for the order parameter. The slopes 
are ^/{2u±) and 7/(2i/||), respectively. The large L limit is easily achieved, obtaining 
7/(2z/x) = 0.96 ± 0.03, in accordance with theory. The other limit is hard to reach, 
because of the slow relaxation in transverse direction. In any case, we found ■y/{2iy\\) = 
0.37 ±0.04. 

The order parameter or susceptibility at has an extra factor which depends on 
the ratio S, e.g., 

*(Tc, L, M) = M-^/^^-^iL^^/^^i /M). (11) 

Moreover, the scaling function obeys ^{S) 5-^11 /(2i^-l) for S ^ 00, and ^{S) ^ 
5'i/2-/3/i'-L foj. 5' — > 0, as a consequence of the limiting behaviors for very long strips 
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[Eq. (9) and (10)] and the hyperscaling relation 



2/? + 7 = zy^ + zy||. (12) 

The scahng form was tested for Ising model. ^-"^^^ Previous applications to anisotropic 
systems were not very successful. Figure 4 is a scaling plot with theoretical values 
of exponents and Tc = 1.41. Similar plots for T = 1.40 and 1.42 show definite deviations 
from scaling. This supports our choice of Tc. The asymptotic slopes for small and large 
scaling variable S = L^/^/M are consistent with the expected value —1/2 and —3/2, 
respectively. 

In Fig. 5 we plot similarly the susceptibility at Tc = 1.41 in scaling form. The 
asymptotic scaling behavior is borne out, 

x(Te, L, M) = M^/^^x{L''-'^\\ /M). (13) 

The scaling function x{S) for large S is x{S) const, and for 5 — > we have x(S') — > 
S^'"^. The data are in full accord with expectations. 

Figure 6 is the scaling plot for the fourth-order cumulant, 

g{n,L,M) = g{L''^/^\\/M). (14) 

Large finite-size corrections are found here. This is also the case for the Ising model. ^^^^ 
These scaling plots are good evidence that at least the exponent ratios P/i'± and 
^\\/^± are in agreement with the field-theoretic values. Scaling off the critical temper- 
ature will determine the exponents themselves. 

5. Anisotropic Scaling away from Critical Temperature 

We expect when T j^Tc, 

*(T, L, M) = L-^/^ii ^ (L^/^ii (T - Tc) /Tc, S) , (15) 

Leung*^^^ proposed a stronger scaling form when S ^ 0. Figure 7 is a scaling plot 
of the order parameter for S — 2~^/^ with Tg = 1.41, assuming the field-theoretical 
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exponents /? = 1/2, i/j^ = 1/2, and v^^ = 3/2. The T > branch (lower part) obeys 
scahng very well. For the T <T(, branch (upper part), deviation from scaling is large. 
It appears that the scaling region for T < is rather narrow. The size effect is further 
complicated by a 1/M correction due to the presence of interfaces below Tg. It is also 
not known how much of this can be attributed to possible logarithmic corrections to 
scaling. Comparing with Leung's data,'-^'^'^-' we feel that his conclusion on data collapse 
is somewhat too optimistic. 

Our 1024 X 64 order parameter data can be compatible with the exponent (3=1/2, 
but only in a rather narrow critical region of AT = 0.05 (see Fig. 8). The data can be 
fitted to a power /3 pa 0.3 in a large temperature region (AT = 0.2). This is also roughly 
the value found in previous work on square geometry. "^^'^^ However, the evidence for 
j3 ~ 0.3 is not convincing since the extrapolated critical temperature would be too low. 

Figure 9 is a scaling plot for the susceptibility, x'? above the critical temperature. 
Good data collapse is obtained with the theoretical exponents. The asymptotic slope 
for large scaling variable M^(T — Tc)/Tc is approximately 7 1. We should not 
worry about the deviations from scaling for large values of T — Tc. Just like the order 
parameter, the data below the critical temperature (not shown here) don't scale well. 

6. Conclusion 

Extensive computer simulation of the driven diffusive model is performed. The new 
data are consistent with the phenomenological finite-size scaling theory with the set of 
exponents derived from field-theoretic model. Thus, this work supports the concept 
of universality of the critical exponents. The order parameter data below the critical 
temperature are somewhat difficult to interpret. The data at and above the critical 
temperature conform to standard finite-size scaling. 
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Figure Captions 

Fig. 1. The locations of the susceptibihty peaks plotted against L~^/^, for rectangular 
systems with S = 2'^^/^ (diamonds, M = 32, 64), S = 1/8 (squares, M = 
16, 32, 64), 5 = 2-8/3 (circles, M = 16, 32, 64), S = 1/4 (triangles, M = 
16, 32), and for the square systems (pluses, M = 32, 64, 128, 256). 

Fig. 2. The fourth-order cumulant g against scaling variable L'^/^(T — Tc)/Tc, here 
Tc = 1.41. The system sizes are 16 x 16 (circles), 128 x 32 (squares), and 
1024 X 64 (triangles). 

Fig. 3. Limiting behavior of the order parameter at Tc- Circles are for L — > oo and 
squares are for M — > oo. The straight lines have slopes 1 and 1/3, respectively. 

Fig. 4. Scaling plot of the order parameter at the critical point, Tc = 1.41. Each set 
of data has a fixed M value, M = 4 (solid circles), M = 8 (open squares), 
M = 16 (diamonds), M = 32 (pluses), M = 64 (open circles), M = 128 
(squares), M = 256 (open triangles), M = 512 (crosses), M = 1024 (up 
triangles). 

Fig. 5. Scaling plot of the susceptibility at the critical point, = 1.41. The symbols 
are the same as in Fig. 4. 

Fig. 6. Scaling plot of the fourth-order cumulant at the critical point. The symbols 
are the same as in Fig. 4. 

Fig. 7. Scaling plot of the order parameter away from critical point. The system sizes 
are 16 x 16 (circles), 128 x 32 (open squares), and 1024 x 64 (triangles). 

Fig. 8. Order parameter to some power, (circles) and ^^^/^ (squares), versus tem- 
perature T. The system size is 1024 x 64. 

Fig. 9. Scaling plot of the susceptibility x' above the critical temperature. The system 
sizes are 8 x 16 (circles), 64 x 32 (open squares), and 512 x 64 (triangles). 
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Fig. 1, J-S Wang 




Fig. 2, J-S Wang 
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Fig. 3, J-S Wang 




Fig. 4, J-S Wang 



10' r 




Fig. 5, J-S Wang 
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Fig. 9, J-S Wang 
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